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ABSTRACT 


An analysis is made of in-plane and out-of-plane vibra- 
tions of a uniform filament which is bent at right angles 
and the ends of which are fixed. Shear deflection and rota- 
tory inertia are neglected. Curve sheets are presented from 
which one may find the first several frequencies of vibration 
in terms of the geometrical and physical variables. An 
important application is to the vibrations of piping systems. 

Also, as a separate project, a vibration program written 
by G. Fink for the CDC 1604 computer and modified by Y. S. 


Kim, has been rewritten for use with the IBM 360/67 computer. 
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A. GENERAL 

in order to assure long time integrity of nuclear piping, 
the applicable codes, ANSI B31.7 -— 1969 Code for Nuclear 
mewer Piping and forthcoming edition of Section ILI of the 
ASME Boiler and Pressure Vessel Code, eal ieeor seismic 
analysis of such piping and FeduUlnrewunal UMe soe Verb y monemGime 
effects of seismically indticed vibration be maintained below 
specified acceptable levels. Generally the problem of 
vibration analysis of the major piping systems is so complex 
Eeies a major computational effort is required and this 
demands digital computer analysis using quite sophisticated 
programs. However, in the early stages of a design, when 
alternate layouts are being investigated, it is useful to be 
able to make some reasonable estimates of what a more elab- 
orate analysis would be likely to reveal. 

A recent document [3] concerned with assisting a designer 
to understand and comply with ANSI B31.7 contains material 
summarizing the vibration characteristics of several simple 
Piping configurations. It suggests how the designer might 
employ this sort of information to assist him in estimating 
the vibration characteristics of more complicated configur- 
ations. This material should prove to be quite useful. 

However, the scope of the configurations about which 


such information is available is quite limited. One of the 








basic configurations with which a designer deals is a simple 
ell-shaped configuration. It was decided that it would be 
quite useful to determine and present in a concise form the 
Significant results relating to the vibration of such a 
configuration. Specifically, it is assumed that the two 
ends of the configuration are perfectly fixed or clamped, 
that the pipe is of uniform properties along its entire 
length and is of rather small diameter compared to its 
memech (hence the word “filament" in the title of this 
thesis) so that shear deflection and rotatory inertia pould 
be neglected, and that the right angle junction of the two 
legs of the ell is small compared to The overall dimensions 
Sueune ell itself, so that it is not necessary to consider 
the properties of the bend or elbow that would exist in 
Bractice at the junction. 

such results are obtained and given herein. The basic 
parameter is the ratio of the length of the shorter leg to 
that of the longer leg; this parameter is designated by the 
Greek letter rho in what follows. A second important 
parameter, Greek beta in what follows, is required to relate 
the coupling that takes place between bending and torsion in 
the case of what are called out-of-plane vibrations. 

This simple ell-shaped configuration has the great 
virtue of being such that the vibration problem is naturally 
divided into two disjoint sub-parts. One such part involves 
only displacements which lie in the plane of the ell; no 


torsion is involved in these cases. The other part involves 








Boch bending andwtorsion, bute hAe mori Onmol (any poimy ome une 
axis of the filament is only perpendicular oe the plane of 
the ell. This gives rise to the two terms used in what 
follows: “in-plane vibrations" and “out-of-plane 


waprations." 


B. ASSUMPTIONS AND LIMITATIONS 

It is well to repeat here specifically the assumptions 

end limitations involved in this investigation. 

iene Conti puraviOneis formed fromeorieginally strarene 
uniform elastic members which meet at a right angle. 

2. The angle at which the two members meet remains a 
right angle throughout the deformation of the 
Conlaeuracgaon., 

3. The ends of the configuration are assumed to be 
perfectly clamped. 

4. The effect of shear deformation has been neglected. 
Only bending and torsional compliances are con- 
sidered and deformations are assumed to be Ss maeeliems 

5. Rotatory inertia has been nomileeead., 

6. Poisson's ratio has been taken to be 0.3. 

Additionally, there are limitations implicit in the fact 

that values of the roots of the equations which are later 
meved are Obtained only to an accuracy of one part in one 
Million parts and in the fact that the Perea eional pro- 
grams used a mantissa having 56 bits. These limitations 


resulted in being able to obtain only about the lowest 





thirteen frequencies for a particular configuration. 


VIL Sore 


only a few values of beta were investigated and there are 


mome difficulties Inet hesuseuol Manear ain berp olay ene t 


other values of beta. 


OG 


A, B, C, D, F,H Constants in equation (2-5) and (2-14) 


q wt Ve Ss FF Oo A 


as 


KS 


WD Q < 


~< 


mah 


NOTATION 


Modulus sot elasticity 


Shear modulus of elasticity 


Functions defined in equation (2-6) 


Moment of inertia of cross section 


IFelLeve omlenny yt Tisiezalel 
Length of pipe 

Bending moment 

Defined in equation (2-4) 
ite 

Torque 

shear force 

Distance from fixed end 
Lateral deflection 

Lateral deflection mode shape 
Slope 

Defined in equation (2-13) 
Defined in equation (2-29) 
Phase angle in equation (2-11) 
Angle of twist 


Angle of twist mode shape 





Constant in equation (2-2) 
Mass per unit length 


Mass moment of inertia per unit length 


Length ratio 


Circular frequency 
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Ti. THEORETICAL DEVELOPERS 


A. BENDING ANALYSIS 
Neglecting the effects of shear deflection and rotatory 
inertia the equation governing lateral vibration of a uni- 


form beam is 


3 4° 
EI rs u=— = 0 (2a) 
OX ot 


fiere Wt 1S mass per unit length, El is flexural rigidity, 
y is lateral deflection and t is time. Assuming harmonic 
Soeillation with circular frequency w, the solution of 


equation (2-1) can be written in the following form: 
Vi ae Xe) COs ict Uy) (2-2) 


mere WY) iS a constant and Y is a function of x alone. sub- 


stitution of equation (2-2) into equation (2-1) yields 


Lt 
ow ~ P'¥(x) = 0 im =) 
X 
where 
2 
ho ow 
ea ET Ba) 


The general solution of equation (2-3) may be written as 
Y(x) = Ag, + BE, + CE. + DE), (2-5) 


where 


dei 








ees eCosh( Px) tucoaGe x (a) 


G5 = %[Sinh(Px) - Sin(Px)] (b ) 

(2-6) 
E25 *[Cosh(Px) - Cos(Px)] (c) | 
g, = %{Sinh(Px) + Sin(Px)] (d) 


At any location x, deflection tL» slope Y_,bending moment 


M and shear force NE can be expressed as 


ix? 
= Ag, + Bg, + CE. + Dey (a) 
ty = P( Ag. + BE 2 ap CE y + Dg,) (b ) 
(2-7) 
7 2 
M, =e bP (Ag. ~ Bey + Ce, + Dg.) Ce) 
- 3 
Ve = EIP (Ag), + Be, + CE, + Dg 3) (da) 
At x = 0, 
Y,=A a. = DP elo 
0 0 (2-8) 
My = C(EIP“) ; Vo = B(EIP>) (c,d) 


Expressing the constants A, B, C and D in terms of YQ» 


ae My and V, and substituting into equation (2-7) yields: 
Y “ EIP3| | y 
Ey g\/P g2/EIP Bay, 
t 2 t 
iy. PE. By &),/EIP E3/EIP Y 
= ; (2-9) 
M Pee E3 EIPg. Ey g),/P M 
V EIP=g, EIP’g, Pg g V 
y 5 2 ii 


x 0 


i 





The above equation expresses the deflection, slope, bending 
moment, and shear force at any location x in terms of those 


at x = OQ. 


B. TORSION ANALYSIS 
The equation of motion of a uniform cylinder undergoing 


mersi1onal vibration is, 
—— = = 0 (2-10) 


where v is mass moment of inertia per unit length, J is 
polar moment of inertia of elastic resistance, and 6 is 
angle of twist. For a sinusoidal solution at circular 


frequency wo, 
CiCxerne— d(x) Com@ne + 7) (2-11) 


where y is a constant and 9 is a function of x alone. Sub- 


stitution in equation (2-10) gives 


O 
qe x) RuRaGe eee (2-12) 
ax 
where 
o 
ac - Vw {Q=) 
JG 


The general solution of equation (2-12) is 
Ox) =F Costex)) tt) sin (ox) (2-14) 


and torque T is 


= soot) = JGa[-F Sin(ax) + H Cos(ax)] (2-15) 
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b) =F 3 TT) = H(JGa) (2-16a,b) 


Expressing constants F and H in terms of 6 and To and 
substituting into equations (2-14) and (2-15), the angle of 
twist >, and torque Te at any location x, can be expressed 


in terms of T, and %% ase to liltows: 


0 
d Cos (ax) Sin(ax)/JGa|i¢ 
It -~JGa Sin(ax) Cos (ax) ib $8 
C. IN-PLANE AND OUT-OF-PLANE VIBRATION ANALYSIS OF AN ELL- 


SHAPED UNIFORM FILAMENT WITH CLAMPED ENDS 


1. In-plane Vibrations 


In Figure 2-1, an ell-shaped configuration and the 


corresponding free body diagrams are shown. 





V 
Moo 02 


Figure 2-l. Ell-shaped Configuration and Free 
Body Diagrams for In-plane Vibrations. 
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The ratio, p; of the lemgth of Che shorverailes tosvear 
or the longer is a basic geometrical parameter in Che 
following analysis. 

At the clamped ends deflections Yoq> Yoo and slopes 
TE ie are equal to zero. Also at the junction of the two 
Bercvions, deflections Try and Yro 
cation of equation (2-9) to the two portions yields, 


id 


are equal to Zero. Appli- 


meer, ea, 0 (a) 
— a 
erp* Sl prp3 Ol 
Bh & 3 t 
Beate ee ae ey (b) 
ey ol eee Lael 
Sy 
em Ol) op Ol, oa Cc) 
(2a) 
53 B. | 
—— M a) = (0 (da) 
gerpS FF prp3 Ye 
By E> t 
—— = i ny = (e) 
Oo Cle ae Ge” “Oe 
: By 
2) op! a Vee eae 2 


where an indicates evaluation of g,(x) for argument (PL) and 
Bs indicates evaluation for argument (pPL). 


At the junction, 


= t e ; = eat 
io [oe i Toone (2-19a,b) 


Combining equations (2-19) with equations (2-18), the 


following set of equations is obtained. 


» 





g3/EIP® 0 g>/EIP> 0 





0 B,/EIP“ 0 é,/EIP> 
_ y 5 =0 (2-20) 
g/EIP -&),/EIP E2/EIP -&2/EIP 


For a nontrivial solution, the determinant, D, of the above 
matrix must vanish. Elementary manipulations on D lead to 


the following form. 


= 0 (2250) 
al 5, =&y By 


The terms of equation (2-21) are functions of p. For 
any value of 9, the zeros (PE), of the determinant can be 
found. The corresponding frequencies are evaluated by using 


equation (2-4) as 


El 
7 (2-22) 


ie = (AL), 
BL 


a 


Solution of equation (2-21) will be explained in 


Section III-A. 


eo. Out-of-plane Vibrations 





Out-of-plane vibrations can be analyzed similarly. 
However in this case the effect of torsion must be taken 


Hmaco account. 
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In Figure 2-2, free body diagrams of an ell-shaped con- 


figuration for out-of-plane vibrations are shown. 





Figure 2-2. Free Body Diagrams of an Ell-shaped 
Configuration for Out-of-plane Vibrations. 


At’ the clamped ends deflections Yo. Yoo and slopes 


NG vl and angles of twist %o1> are equal to zero. 


t t 
mi *02 %O2 
Application of equation (2-9) to the two portions yields 


IL 





g g 
Ye Ey (a) 








een ee 

Ya = = “ony a tb) 

tg Si eak “ Yon Ce) 

Vial Seoon Baga (d) 

oo Bs i B. : al (2-23) 
EE MODAL ae eNO 

eee = Mast 1, liog (f) 

Mio = Moo + oe (g) 

Vio = P&oMos + 81Vo0 (h) 


where g, indicates evaluation of g, (x) for argument (PL) and 
gE. indicates evaluation for argument (pPL). 


For the torsion, equation (2-17) applies and gives, 


_ Stal eil) 
Oe ae ee (a) 
fl JGa Ol 
Tr, = Cos (aL) To Go, 
(2225) 
_ Sin(apL) 
nD = TE 62 Cc) 
Tro = Cos(apL)T,. | (ad) 


por Mem iuMmecus Ol: 
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M = T : T = —M (a,b) 


t1 7 The mI L2 

intel : ot _ 
ee ee Pee Ss (c,d) (2-25) 
Voie Mp 8 a Ee (e,f) 


specializing to the case of circular cross section and 


ecsuming Poisson's ratio is 0.3, then 


—— = = (2-26) 


Combination of equations (2-23), (2-24), (2-25) and 


(2-26) givesthe following set of equations. 





BP “E3P 0 me Bo 
2 1, 3P¢ 
0 EP - s—Sin( OL) 0 0 
P 0 0 1.3?“ sin (apb) 
Ey 0 63 
— aa. 
eee BP 0 q oa 
Sy 
Ey 0 0 -Cos(apL) = 
0 g, Cos (aly 0 0 
( 22278) 


For a nontrivial solution, the determinant D of this matrix 
must vanish. From equations (2-4), (2-13) and (2-26) it is 


meund that 


(aL) = B(PL)* (2-28) 


1, 





where 





B A 5 (2-29) 


is a dimensionless parameter. Substituting equation (2-28) 
into equation (2-27) and performing elementary manipulations 


on D give the following equation 


B36 : : Bo 82 

0 By ->52Sin[6(PL)*] 0 0 Bs 
0 0 =sSsin[pB(PL)“] gz 0 |= 0 
g = = 

2 5 0 0 By ay 

(2-30) 

0 0 (PL)Cos(Be(PL)"] gy 0 

Me &, (Pi)Cos(e(PL)“ 0 re, 


The terms of equation (2-30) are functions of p and 8B. 
For given values of pep and 8, the zeros (PL), of determinant 
D can be found. The corresponding EreauencToemercc evaluated 
by using equation (2-22). 

Solution of equation (2-30) will be explained in 


section ILI-A. 
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Lit > se PRE OE NGA IO NSC ieee Ue 


fm METHODS OF SOLUTION 

The equations (2-21) and (2-30), developed in Chapter II, 
can be solved either by algebraic manipulations or by using 
Girect computer solution. 

In a first approach to a solution, algebraic manipula- 
tions were performed to reduce the order of determinants. 
Then transcendental equations were obtained and two computer 
programs, one for the in-plane solution and one for the out- 
of-plane solution, were developed to get the roots of these 
equations. 

In the computer programs the transcendental equations 
were used as function subprograms and roots located in an 
interval. Then the Newton-Raphson method was used to obtain 
good approximate values of the roots. 

The procedure explained above led to some difficulties. 
During the-solution of the equation which was derived from 
equation (2-30), the Newton-Raphson method failed to give an 
expected root because the slope of the transcendental function 
reached very high values and this caused a jump to adjacent 
roots. In one case, the number of iterations required 


8 (for 


6 


exceeded three hundred with an allowable error of On 
the quantity PL) and increasing the allowable error to 10 
did not make much difference. Finally, comparison of 


results with the classical solution in the case o is 


ail 








equal to zero showed that some of the roots were 
exc raneous. 

The difficulties mentioned in the preceding paragraph 
were eliminated by making direct computer evaluation of 
the determinants (2-21) and (2-30). For this purpose 
Subroutine BRODE was developed and each element of the 
@eterminants was stored in this subprogram. In addition to 
emprouctine BRODE, library subroutine DTERM which was 
developed by Jean Bow in 1967 was used to evaluate the 
determinants. This modified procedure resulted in reduction 
of computer time required and in the elimination of extra- 
meous roots. 

Perr CultTes With OCcastonaleiarlume of the Newvon— 
Raphson method were reduced by using the method of successive 
bisectioning. In this method the number of iterations 
required for one root did not exceed fifty with an allowable 
error of Ono This procedure also resulted in conservation 
of computer time. 

Detailed instructions for use of the computer programs 


are given in Appendix B. 


B. GRAPHS 

Numerical results given to a large number of significant 
figures were obtained primarily for the purpose of investi- 
meine the aceuracy and integrity of the solutions. Chapter 
IV will give a discussion of these investigations. 

However, so as to present the results in the most con- 


femreny Lorm for practical use, graphical outputs 


Bia 





were also obtained by use of the library subroutine 
DRAW. | 

One graph was obtained for the in-plane solution. For 
the out-of-plane solution five graphs were obtained cor- 
responding to five different values of 8. In Senn Cecee 
the ordinates are values of (PL) and abscissas are values 
of op. On these graphs the first (lowest) curve is for the 
first mode, the second curve for the second mode and so 


on. Graphs are given in Figures 3-1 through 3-6. 


C. SAMPLE PROBLEM 
A sample problem will be solved to show how to use the 
Praphs Given in the preceding article. 


Properties of ell-shaped configuration: 


Ly kG’ —- Te 

Li LO mn oe 

Outside diameter of pipe : 6.625" 

Inside diameter of pipe eo OOS. 

Wall thickness of pipe ae OPN OTe 

Weight of pipe a ule lem Ihey Ais 

Micekness of insulation ye 

Weight of insulation ee OOD fae 

Weight of contents 5 LZ la loy ey 

Modulus of elasticity E : 27,400,000 1b/in® 
G : 10,600,000 1b/in“ 


For the use of graphs, values of o and 8 must be known. 


mor this problem, 
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Values of Rho 


In-plane Vibrations. 


24 





Values of (PL) 






| 
0.4 


Values of Rho 


Figure 3-2. Out-of-plane Vibrations, Beta = 0.00001. 





me inns 
See eciie 
CAO 
PAT 
B2ernini 


oe) oe) MN 








(Id) go sente, 


Values of Rho 


Out-of-plane Vibrations, OE O05. 


Figure 3-4. 


aad 





Values of (PL) 


Bere 3-5. 


Values of Rho 


Out-of-plane Vibrations, 
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8 = 


o- 





Jake 





ile) 


ULE 
HU | 


0.8 





(Td) JO SoOnpTen 


Values of Rho 


Out-of-plane Vibrations, CreC 5: 


Figure 3-6. 


aa 








je /esmied Cam 


where v is the total mass moment of inertia per unit length 
for the torsion and uw is the total mass per unit length 
undergoing lateral vibrations. Assuming that the contents 


(water) do not rotate, then 


, + ve : 
Ypipe Vinsulation 


tov & Moipe vi Vinsulation i Yeontents 
numerical values of Yuin and We otal are 
= 2 
aa 0.0509089405 lb-seec 
= Oo. 2 
en 1.068875490 l1b-sec”-in 


and from equation (3-1) 


B = 0.0150806896 


For evaluation of frequencies, equation (2-22) should be 


used as 
a 2 aa 
W. = ae EL 
pL 
3 
then 
= = Qe Den Sme Gi eee 
HL, 


30, 





To solve the problem, values of (PL), can be found from 
the graphs. But no graph is available for the value of 8B 
found above. Therefore interpolation is required between 
the values taken from Figure 3-5 and from Figure 3-6. Values 
of (PL), are given in the table below where the exact 
solution was obtained by a computer solution using the 


specified value for B. 


MODE EXACT SOLUTION OBTAINED BY 
NO. SOLUTION LINEAR INTERPOLATION | 
il 2.624579 2.624302 
2 isos 4.521991 
3 6.729044 6.706905 
4 7.862686 i. Oleigee 
5 TO) g Byles Gs Nels Bsior 22 
6 11.9087 33 11.698410 
7 N35 OP! Ty 13.442683 
8 14.259158 15.619057 
9 16. 376500 16.379861 
10 17.185014 17. 403866 


As seen from the table, for the lower modes, linear inter- 
polation gives good results. However there is greater error in using 
linear interpolation for the higher frequencies. One can 
make what appears to be the reasonable conjecture that if 
one were to plot curves with (PL) as ordinate and 8B as 


abscissa and with different curve sheets corresponding to 


oul 





fixed values of pop, there might be some crossing of individual 
curves; that is, for certain values of B Somatic modes 
might have the same frequency. Clearly this would confuse 
The significance of characterizing different modes simply 
mamordinal numberwain thewerder of inereéasrinee requcncy sand 
mould indicate that a more sophisticated interpolation is 
meomared.,. However, no such graphs were prepared and there 
do not seem to be enough values of 8 represented in the 
graphs at hand to permit making replots. Thus it is 
suggested that some future investigator might find it 
interesting to employ the computer programs, incorporating 
eect modifications, to obtain such graphs and thus inves-— 
tigate the reasons for inaccurate interpolation for the 


higher modes. 


OU)’ 
N 





IV. ACCURACY AND ASSURANCE OF SOLUTION 


The allowable error permitted in determining the roots 


6 


M10 and this was assumed to be satisfactory for practical 
purposes. The allowable error can be changed from 107° ZO 
any desired value by changing the value EPS in the computer 
programs. This can be done by changing the card 000090 intre 
program IN-PLANE and card 000100 in the armenian OUT-OF-PLANE. 
Assurance of solution was established in the following ways. 
1. In the case p is equal to zero, the ell~shaped con- 


figuration becomes a simple built-in beam and solution for 


this case can be compared with the classical solution. 


Degenerate Case; 


Vices. Ou nt) 


Comparison In - plane Out-of-plane 
Values * Soul Uipal eye Sie sLphasl@jla 


ul 
2 
5 
U 
5 
6 
i 
8 
? 


[oJ 
2 


% 
Wbtained by computer solution 


4.7300407 
7.8532046 
9956078 
a LS ess 
-2787597 
PAO 5 22 
5619449 
- 7035376 
8451302 
.9867229 


mesa (PL) Cos(PL) = 1. 


7300289 
8531849 
9955803 
micyaeoe 
2 {eyes 
4203012 
.5618858 
(rOsto14 
nod (ae 
9890104 


of the equation 


4.730041] 
foes ae iy 
9956081 
HAS) (esis 
.2787594 
4203526 
5619449 
5 OS. 
.8451861 
.9866973 








In equation (2-21) when p is set to zero the value of 
the determinant is zero regardless of the value of (PL). 
Therefore p was taken to be 0.00001. 

in out-of-plane case for the larger values of 8, pure 
torsional modes appeared and this can be seen from Figures 
(3-5) and (3-6). The values of (PL) given in the above 
Cable were taken when 8 is equal to 0.001. 


2. Consider the cases shown in Figure 4-1. 





(a) (b ) 


Figure 4-1. Cases Having Reciprocal Values of op. 


The configurations shown in the above figure are identical, 
therefore two systems must have the same frequencies. Using 


equation (2-22) as 


(ame (Er (na (EX 
ou rea § 


then 
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For out-of-plane solution in the case 8 = 0.001 


No. Co (PL) 
bi 


2.9832348 
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3 Be OAC IE 3.6886158 -999993485 
4 .8416393 4208958 -999965571 
5 .8159008 5. 4080044 .999980037 
6 6035633 6.8018597 -999977062 
7 .0968311 7.5488148 .999 894233 
8 . 0625606 8.5315098 .999946210 
S . 8673419 9.9339245 -999948965 
10 » 3568783 .6796268 oo oe 
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3. In the case p is equal to one, it is expected that 
among the in-plane solutions will be many for which the 
moment at the joint vanishes. Thus, these cases also per- 
tain to the clamped-pinned simple beam. Now two clamped- 
pinned simple beams of different lengths, &£ and 2*, may have 
Mevching frequencies in their spectra. That is, for example, 
the second frequency for one may equal the first frequency 
memecne other. Specifically the following case was 
wensidered. 

The third frequency for an equal-leg ell is equal to 
the second frequency of a clamped-pinned beam of length &, 


namely 


2 EL 
Cr 
3 Say) nt 
where 


(PL) 3 = 7.0608586275 


The first frequency of an equal-leg ell of length 2* is 
equal to the first frequency of clamped-pinned beam of 


length &£*, namely 


_ 2 EI 
of Niles ELL 


where 


(PL), = 3.92660231 


The ratio P was determined so that then the frequencies are 


equal, thus 


BE 





Then performing an analysis for an ell-shaped configuration 


with the value of op, found above, the second frequency is 


indeed, 
Ms Dal 
Wo = (PL), \ [7 
ys 
where 


(PL). = 7.068384686 
This value shows close agreement with the value 
(PL). = 7.0685275 


given above. 
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APPENDIX A 
DESCRIPTION OF PROGRAMS 


A. GENERAL REMARKS 

Programs IN-PLANE and OUT-OF-PLANE are FORTRAN-IV 
language digital computer programs designed for use in the 
vibration analysis of an ell-shaped uniform filament with 
clamped ends. 

Both programs employ double precision arithmetic. Sub- 
memeaines are provided in binary form. Trigometric and 
hyperbolic functions, available in the IBM-360/67 library, 


were used. 


B. PROGRAM STRUCTURE 

Each program consists of a four section main program and 
tree subroutines. The function of each SCCUlCN OL oe em arna 
program and of each subroutine is given below. 


1. Main DreOcranmeupseoculonc and  —vnel rar uUnet Tons. 


ENP UT Controls read-in of data and allocates 
SvOrace slOcCauLOnseiOlr aterays and SUubSCripred 
Vieiel apes es Nera! Keane “error Tor Gever= 
Mination of roots is specified in this 
Section. 


mle RATION In this section the roots of the determinant 


AND are located in an interval by searching the 
EVALUATION sign changes of the determinant and then the 
OF ROOTS method of successive bisectioning is used to 


evaluate the roots. The maximum number of 
iterations was limited to three hundred and 
this was assumed sufficient. If the number 
of iterations becomes as large as three 
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CErPur 


GRAPH 


hundred the process terminates and the 
approximate root is tagged to indicate a 
possible inaccurate result. 


Controls output format. 
Controls the graph output and uses the DRAW ~ 


subroutine provided in the IBM-360/67 
iporary . 


rae SUDPOUTANeES and thelr Tunceviensc. 


BRODE 


DTERM 


In program IN-PLANE equation (2-21) is stored 
and in program OUT-OF-PLANE equation (2-30) 
is stored. Evaluates each element of the 
determinant and obtains the value of the 
determinant by calling subroutine DTERM. 


This subroutine computes the determinant of 
a matrix of real numbers by Gauss‘! method 
of elimination with row pivoting. 


PROGRAM NOMENCLATURE 


ne), nA 


Elements of determinants. 

Checks sign of determinant and if there is a 
sign change moves control to the method of 
successive bisectioning. 

Values of 8 in Subroutine BRODE. 

Checks the acceptability of roots. 

Veleibisy ene ts) 

Single precision version of RO. 

Value of determinant. 

Allowable error in roots. 

Value of determinant. 

Index of modes (Ordinal number). 


Values of (PL) used to search interval. 


Valuccectedercecrmanane inthe Metnou OL 
successive bisectioning. 


Values of 0. 





ROOT 


SEV 


Ros 


1a 


XX2 


Roots of determinant (PL). 
Corresponds to RO in subroutine BRODE. 
Sulighelle sjonceres Si lieial Wdcietsialelel end IsOlOML 


Middle of interval. 


Upper and lower limits of interval contain- 
ing the root. 


Value Ss Or eP awk. Clondeoe2 in sb rougaie 
BRODE . 
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APPENDIX B 
INSTRUCTIONS FOR USE OF PROGRAMS 


A. DATA REQUIRED 

ieee Lata reguaredsicr prepram) IN-PLANE 
a. Number of o's (NRO). 
b. Number of modes sought (MODE). 
c. Graph control (NGRAPH). 
d. Values of 9 (RO). 
é.. Graph title (LTMBEE)e 

meee Data reguired fommerogram OUT-OFP-FLANE 
a. Number of o's (NRO) 
b. Number of 8's (NBETA). 
ec. Number of modes sought (MODE). 
ad. Graph control (NGRAPH). 
2. Welles wor (2 (islaZbe 
f. Values of p UROL) 


ge. Graph title (ITITLE). 


ee DATA DECK ASSEMBLY 
1. Program IN-PLANE 
The data deck is ordered as follows: 
a. Problem statement card. 
b. Values of so cardcs). 
CG. - Graph Gitlescards: 


2. Program OUT-OF-PLANE 





Whe data deck 2s .ordered aoc lovs. 
a. Problem statement card. 

b. Values of B card(s). 

c. Values of p card(s). 


ad. Graph title cards. 


C. LIMITING THE NUMBER OF ITERATIONS 

In the programs the maximum number of iterations was set 
fen500. If the number of iterations reaches 300 for any 
root, the current value is assumed as a solution but in the 
output it is marked as ITER. EXCEEDED 300. A user either 
accepts the root or changes the cards 000430 in program 
IN-PLANE and 000460 in program OUT-—OF-PLANE by changing 300 


to any desired value. 


eee GRAPH 

A user can get graphical output by option. DRAW sub- 
routine, available in IBM-360/67 library was used and three 
CALL DRAW statements control first, intermediate and last 
curves to be plotted on one graph. DRAW subroutine requires 
single precision arithmetic. Therefore variables CEM and 
SEV are introduced to change the values of RO and ROOT from 


double precision to single precision. 


BR. FORMAT OF DATA 


1. Program IN-PLANE 
a. Problem statement card. 
The card has 3 fields under the control of 


following field specification. 


HO 





313 
Field names and the information each field conveys are 
Biven below in order. 

Field names and content. 

(1) NRO Number of o's 

(2) SNMODE Number of modes sought. 

(3) NGRAPH Code specifying whether there is 

Co be srankheeulpur- 

(NGRAPH = 0 No graph output) 

(NGRAPH = 1 Graph will be output) 
b. Values of o card(s). 

Hield speciirtearion (HEC s0" 

If there are more than 7 different values of 
op, additional card(s) will be required with the same field 
specification. 

Coreareaph wat le cards. 

The title of graph is read as data with the 
field specification 6A8. For the title two cards are 
Mmeguired,. ithe first card corresponds to the first line of 
the title and the second card corresponds to the second 
line. 

2. Program OUT-OF-PLANE 
2, (eile lam Breelieeileinn Genes 

The card has 4 fields under the following field 

specification 
4J 3 


Field names are given.below in order. 
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(1) NRO Number of o's. 

(2) SNBEL Number of Bis, 

(3) MODE Number of modes sought. 

(4) NGRAPH Same as in program IN-PLANE 

Wellocles, (ce (6) (eels 

Field specification 7F10.0 and field name is BETA. 
Values of 09 card(s). 

Same as in program IN-PLANE. 

Grapn Citta cards. 


Same as in program IN-PLANE. 


4 


ALPEN DIAC 
FLOW DIAGRAMS 


A. GENERAL 
Flow diagram of the program OUT-OF-PLANE is given. This 
flow diagram includes all parts of the main program. Flow 
diagrams of subroutines BRODE and DTERM are not given and 
their functions were expldined in Section A-2 in Appendix A. 
Program IN-PLANE is similar to OUT-OF-—PLANE, therefore 


its flow diagram is not given. 


B. FLOW DIAGRAM OF PROGRAM OUT-OF-~PLANE 


Flow diagram follows. 


ae 








READ: NRO,NBETA 
MODE, NGRAPH 









READ: BETA(II) 
hee een 






READ: RO(JJ), 
Jd) = Re 


Ii 







CALL BRODE(P(2), 
REG). 
BeAr) bry) 





CALL BRODE(P(K) 
Qi! ji - 


BE aCe ee DE) 
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He) es DIRE 


B=F(K)*F(K-1) 
< 

Pa 

ia 


a0 






Xa 72 
+(XX2)/2 


BBB = XX2 -— XXXl 


< 
_— ROOT(N,J,L 
BBB:EPS = eee 


> 


CALL BRODE 
COL HOG je 


BETA(1),DET) 





R3 = BET 


CALL BRODE 
(X(I) ,RO(T), 


BETA(L) DET) 





RY = DET 





ne 
P(K+1) = 





> 
CS i I <:300 > 


WRITE: ROOT 
(N,J,L), ITER. 
EXCEEDED 300 





ae 


P(K) + 0.05 


J+1 <i 
ea 


WRITE: 5 
ROO Mma one) 
1 ,MODE 


1,NRD 
1, NBETA 


NGRAPH:0 > 8 


a 
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READ: ITITLE(T) 
faa=) eG 
READ: ITITLE(T) 
ie ae 







MM = 1 
MMM = 1 











en = 
ROOT (M,MMM ,MM) 


<a CALL DRAW 


CALL DRAW 


CALL DRAW 


NM: NBETA 





AS 





<itinsert fn sae 
a. 
<finone > [a+ 


SLU 


END 
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APPENDIX D 


PROGRAM LISTING 


A. INDEX 


1e8 


ae 


Program IN-PLANE . 
Main program: 
section 
INPUT 
ITERATION AND 
EVALUATION 
OF ROOTS 
OlORL Ie Uk 


GRAPH 


Subroutines: 
BRO DE 


DTERM 


Program OUT-OF-PLANE 
Main program: 
section 
INPUT 
ITERATION AND 
EVALUATION 
OF ROOTS 
OUR ea, 


GRAPH 


Sequence number 
000090 —- 000130 


000170 - 000550 


O005705— 0007 70 


000810 - 001030 


001060 — 001300 


001310 - 001610 


pequemcer WUIlber 
000100 - O0O0L50 


000190 - 000590 


000630 - 000830 


000870 - 001110 


Sl 


5D 
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Subroutines: 


BRODE 


DTERM 


LISTING 


001140 — 001600 


001610 - 001910 


ae 


D7) 
60 
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PROGRAM QUT-OF-PLANE 
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APPENDIX E 
PROGRAM VIPIPE 


A. GENERAL 

As a part of the present thesis investigation, attention 
was given to a digital computer program called VIPIPE which 
was originally written by Fink [1] and later modified and 
augmented by Kim [2]. This program provides solutions for 
in-plane and out-of-plane vibration of a single plane piping 
configuration having quite general geometry (except, of 
course, for its limitation to a single plane). It was 
desired to convert this program from a FORTRAN 63 program 
for use with CDC-1604 computer formerly employed at the 
Naval Postgraduate School, to a FORTRAN IV program for the 
IBM 360/67 presently in use, in order to provide an indepen- 
dent means of checking the results obtained as indicated in 
the body of this thesis and, also, because it would be 
useful to have a program of such generality available. Por 
details of the theory of VIPIPE, refer to the theses written 
by Fink and Kim. 

Because of the limitations on word length in IBM 360/67, 
the converted program, which was called VIPIPE-1, has been 
Veerueen tO employ double precision arithmetic. Certain 
improvements have been made in the capability to produce 


prcaphical output. 
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The effort to make this conversion has been found to be 
much greater than was anticipated. It seems that CDC-1604 
computer is a more "forgiving" machine than is the 
IBM 360/67. In checking out the new program certain bugs 
fame FO light which were actually present imetiie eceriier 
version but which apparently did not cause trouble. However, 
they represented some syntactical errors that resulted in 
failure on the IBM 360/67. In some cases it was clear how 
tO modify the program so as to eliminate the errors but in 
other cases the difficulty was so deep in the logic that a 
Serreccion could not be made without acquiring as deep a 
knowledge of the program as was possessed by original 
Mereers, Fink and Kim. 

Thus it cannot be asserted that VIPIPE-l has been 
successfully converted in every fine detail and it is 
recommended that additional efforts be made to assure the 
overall integrity of the program and to assess its limita- 
tions. It should be noted, however, that no gross errors 
have been found in any results obtained by use of the new 
meeemam; the only question is that of the degree of accuracy 
taae it is possible to obtain. 

femparisons of accuracy are difficult to make. In the 
first place, the word-length is different in the new progran, 
and, other things being equal, one would expect the new 
mMeecran tO provide better accuracy, and to maintain its 
integrity in the face of round-coff error beyond the point 


where the earlier programs "bombed out." Thus, agreement 
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between Kim and Fink and disagreement with VIPIPE-1 does not 
mean that the former are "correct" and the latter is in 
error. Furthermore, Kim reports having discovered and cor- 
rected some bugs which he found in Fink's version. Also, 
there is some difficulty in understanding the precise 
details of some of the problems treated by Fink and Kim. 
With this preface which admonishes the reader not to 
make superficial inferences from comparisons, in Section C 
some comparisons between results obtained by Kim and results 


obtained using VIPIPE are presented. 


B. DATA DECK ASSEMBLY 
The order of data deck is 
1. Graph title data cards. 
25. lemoo lel Sweweme me weuetel- 
3-5 Bolladary conditrven data card(s) - 
4. Branch joining angle data card(s). 
5. Component extensive properties data card(s). 
6. Component intensive properties data card(s). 
ine lOespecinteabtomrand formageo: —2rapn vivtie ere are 


given in Section E-l-c of Appendix B and information for the 


remaining cards are given in [1]. 


Cee OAMPLE PROBLEMS 
Properties of piping systems analyzed in this section 
are given in [1]. Here, only numerical values of frequencies 


will be given. 
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1. System 1 (in-plane) 

Mode Comparison* VIPIPE** ae ir ela 
No Values Delta m. (D) Delta m. 

oil 75.09994540 75.09994544 75.09994535 
2 207.01589137 207.01589036 207.01589099 
3 aOS = SSSI) Ley 405.8339 3183 405.83391910 
4 670. 86408383 670. 86310856 670.86408365 
5 1002.15520225 1002.19978814 1002.15520252 
6 1399.70436269 1399.70434917 
fe 1863.51172599 1863.51198134 


*Sources of comparison frequencies for system 
are given in Section D. 


**Frequencies are given in [1]. 


**¥* Frequencies are found by VIPIPE-1l. 


ZG 


Mode 
No. 


NO 


—[-_ WW 


“NI OW UI 


System 2 (in-plane) 


Values 


Abe Clove S15 ete 
WS oe DO21 e228" 
207.09776595 
405.82896584 
670.86435913 
1OO2 5 15 Le 70 


1399. 70436229 


Viet es 
DVS Te 2) jl, 1D) 


11. 80213389 
73.96272167 
207.09776685 
405.82894132 
SHO ge 2 Olers i 
1002.26833752 
1399.088495 34 


65 


ances ys pele c 


VIPIPE-1 
Delta m. 


Comparison 


11. 80213432 
(Sj 2V 2081 
207.09776559 
405.82896568 
670. 86435889 
NOC, 2 5 Sys auton fis 
1399. 70437163 





3. System 4 (in-plane) 


[a LL 


Mode 
No. 


1 
2} 


3 


ViRUPE 
M-D m. (D) 


66.69934902 
171. 22352538 
S355) 6 SOO Cos 
535.64529414 
773.44358487 


4, System 5 (in-plane) 


VIPIPR=2 
M—-D m. 


67.31642169 
172.86857155 
33909491357 
544,08929568 
794 57343086 


Mode WIL Ie WOES eee 
No. Delta m. (D) Delta m. 
1 116. 89479180 118.95185774 
2 304.88354398 B02 .912005 56 
3 566.43670976 519.19656352 

Mode VIPIPE VIPIPE-1 
No Pave ec) OD) PM ine tb 
1 Mish Lise eis 
2 302.91288592 

519.19656389 

\ 826. 88281929 794.49003694 
5 1126. 83797508 1126 27858174 
6 USNS 4 Syoteny ols |Z 1464.59530724 
i 2099528440 30 1874.97923215 
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D. SOURCES OF COMPARISON FREQUENCIES FOR SYSTEM 1 AND 
SYSTEM 2 


The governing fourth order equation of a homogeneous 


Single component straight system 


LV noe 


ES ae (ES) 


was solved to give 


VecmeaCOsS! X -GmbmmCOSl) f unmiOm 5) Kimenex ee meGier > 1 1) iteme om 


Upon applying boundary conditions for system 1* Grey = 0 
at x = 0, x = L) and simplifying, it was found that the 
eigenvalues of the system must satisfy 


Cos(PE Cosht(PlL) = 1 (E-2) 


Upon applying the boundary conditions for system 2 (y=y' = 0 
at x = 0, y" = y™= 0 at x = L), it was found that the 


eigenvalues of the system must satisfy 
Cos(PL) Cosh(PL) = -1 (E-3) 


The roots of equation E-2 and E-3, from which the natural 
frequencies of corresponding systems may easily be found, 


were obtained through iteration using a digital computer. 


*System 1 is clamped-clamped end beam 
System 2 is clamped-free end beam 


67 





ee FPROGRAM LISTING 
1. Index 


Main program: 


68 


SSG EALole sequence number ene. 
INPUT 000010 - 001080 70 
INVARIANTS 001120 - 001520 72 
CONTROL IN PLANE 001560 - 002470 (ms 
CONTROL OUT OF PLANE 002520 - 003420 5 
ITERATION 003460 - 004460 77 
OUTPUT 004490 - 008135 79 
Subroutines: 
SUBSEC 008140 - 008640 87 
DISTM 008670 - 009220 88 
DISTMO 009652 - 009676 89 
RIGID 009652 — 009676 eee 
RIGIO 009679 - 009730 gk 
STIFCO 009760 - 009970 92 
STIFOO 010000 - 010230 Je 
STAVEC 010260 - 010550 ce 
STAVEO 010580 - 010830 ge 
INVERT 010860 - 011570 94 
HANGER ISO = OUiee 96 
HANGEO TsO = Culeyide 96 
CFIELO 012000 - 012520 oe 
POINO Glass = Clavie 98 
CFIELD 012000 - 012520 98 
POINT 012550 - 013490 Ne) 





Ze 


SFIELD 


SFIELO 


MATMUL 


FINBRA 


FINMAT 


BRANCH 


BRANCO 


STATEM 


STATEB 


BRCOR 


DELMA 


Listing 
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0135120 
Oars OO 
014210 
014360 
014650 
014950 
015610 
016320 
016870 
017630 
0 


013880 
014180 
014330 
014620 
014920 
015580 
016290 
016840 
017610 
017940 
019380 


Oo 
ILOM 
menE 
102 
Oz 
103 
104 
106 
107 
108 
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PROGRAM VIPIPE-1 
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